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ABSTRACT 

We model the annihilation of relic positrons produced in a gamma-ray burst (GRB) after its afterglow 
has faded. We find that the annihilation signal from at least one GRB remnant in the Milky Way galaxy 
should be observable with future space missions such as INTEGRAL and EXIST, provided that the 
gas surrounding the GRB source has the typical density of the interstellar medium, < 1 cm -3 . Three 
fortunate circumstances conspire to make the signal observable. First, unlike positrons in a standard 
supernova, the GRB positrons initially travel at a relativistic speed and remain ahead of any non- 
relativistic ejecta until the ejecta become rarefied and the annihilation time becomes long. Second, the 
GRB remnant remains sufficiently hot (T ;> 5 x 10 5 K) for a strong annihilation line to form without 
significant smearing by three-photon decay of positronium. Third, the annihilation signal persists over 
a time longer than the average period between GRB events in the Milky Way galaxy. 

Subject headings: gamma rays: bursts - supernova remnants - X-rays: ISM 



1. INTRODUCTION 

The detailed nature of the central engine of gamma- 
ray bursts (GRBs) is still unknown (van Paradijs, Kouve- 
liotou, & Wijers 2000; Piran 2000, and references therein). 
Attempts to derive empirical constraints on the environ- 
ment (Djorgovski et al. 2001, and references therein), the 
frequency (Schmidt 2001), the collimation angles and en- 
ergy output (Frail et al. 2001; Freedman & Waxman 2001; 
Panaitescu & Kumar 2001b; Piran et al. 2001) or the pos- 
sible association of GRBs with supernovae (Bloom et al. 
1999; Kulkarni et al. 2000; Reichart 2001) are compro- 
mised by the difficulty of observing GRBs across large cos- 
mological distances. Obviously, identification of old GRB 
remnants in the local universe would provide much better 
insight into the nature of GRB progenitors (Loeb & Perna 
1998; Woods & Loeb 1999; Perna, Raymond & Loeb 2000; 
Paczynski 2001; Ayal & Piran 2001). 

Given the recently inferred similarity between the en- 
ergy output in supernovae and GRBs, ~ 10 51 ergs (Frail et 
al. 2001), it now appears difficult to separate the late evo- 
lution stages of their remnants hydrodynamically. Even if 
a GRB explosion is initially highly beamed, the asymme- 
try of the blast wave it generates in the surrounding inter- 
stellar medium (ISM) would be erased on an isotropiza- 
tion timescale ^ so ~6x 10 3 (i?5i/7io) 1 / 3 yr (Ayal & Piran 
2001), where E§i is the kinetic energy output of the GRB 
in units of 10 51 erg and uq — pism /Vt> TO p is the ambi- 
ent number density of atoms in units of cm~ 3 for a mean 
molecular weight = 1.4. Subsequently, the blast wave 
expands just as in a normal supernova (SN) remnant. Any 
non-relativistic, SN-like ejecta would moderate the initial 
GRB asymmetry as soon as it overtakes the decelerating 
GRB shock (Piran & Ayal 2002). External anisotropy 
may also result from the interaction of the remnant blast 
wave with a non-uniform ISM. Thus, morphological stud- 
ies alone cannot unambiguously identify GRB remnants. 
The alternative method of seeking spectral signatures of 
the photoionized regions around GRBs (Perna, Raymond, 
& Loeb 2000) also suffers from potential confusion with 



SNe occurring in unusual environments. 

In this Letter, we point out that positrons produced dur- 
ing the early relativistic GRB phase offer a powerful tool 
for identifying GRB remnants. We model the expansion of 
a GRB remnant and show that a well-defined annihilation 
line should be visible until radiative cooling becomes im- 
portant, at which time the remaining positrons annihilate 
rapidly (see also Dermer & Bottcher 2000). 

2. MODEL OF GRB REMNANTS 

We model the initial GRB explosion as two highly- 
relativistic jets of material moving in opposite directions, 
initially covering a fraction /& of a sphere surrounding 
the GRB. We assume that a fraction £ + of the GRB 
energy is transformed into e + rest mass. The positrons 
can either be produced during the burst itself (Cavallo & 
Rees 1978; Shemi & Piran 1990) or by interactions of the 
7-rays with the ambient medium (Thompson & Madau 
2000; Dermer & Bottcher 2000; Meszaros, Ramirez-Ruiz, 
& Rees 2001). We scale our results to a fiducial value 
of £_2.3 = (C+/0.005) corresponding to half of the GRB 
energy in positrons with a bulk Lorentz factor 7 ~ 100 
as required by the compactness argument (Piran 2000). 
We therefore assume a total of N + = 6 x 10 54 £_2.3-E r 5i 
positrons in the GRB. 

The positrons initially travel at a relativistic speed just 
ahead of the GRB jet. The jet begins to decelerate when 
it reaches a radius Rj ~ 0.2(i? 5 i/no) 1 / 3 pc (Rhoads 1997). 
At this point, it begins to expand sideways and slows ex- 
ponentially to non-relativistic speeds. Thereafter, sim- 
ulations have shown that Rj oc i 1 / 3 until the remnant 
isotropizes at ij so and begins a Sedov- Taylor spherical ex- 
pansion (Ayal & Piran 2001). 

If the commonly hypothesized link between GRBs and 
SNe exists, then the jet will be followed by SN ejecta 



of mass M P 



10 Mq. The ejecta expand freely be- 



hind the jet and eventually impact the (decelerating) 

positrons at a time t m - 1K ~ 110(i?5i/7io) 1 ^ 3 w 4 3 ^ 2 yr, where 
«4 is the ejecta velocity in units of 10 4 km s _1 (Piran 
& Ayal 2002). At this point, the density is n m ; x ;> 
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50Afi(rio / E$i)v 4 cm -3 , assuming that the ejecta are dis- 
tributed uniformly throughout the remnant and where 
M\ = (Mej/10 Mq). If mixing of the ejecta and shocked 
gas is efficient, then the ejecta-dominated expansion mod- 
els of Truelove & McKee (1999) show that by this point the 
thermal energy imparted to the swept-up material could 
heat the remnant to a temperature r m ; x ^ 10 8 K (E 5 i/Mi). 
The ratio of the annihilation time to remnant age does not 

change with time, r ann /i ~ bOOM^ 1 {E$i / n a )v± Z ^ 2 . Thus, 
only a small fraction of the positrons should annihilate 
during this very early stage, even if the ejecta are concen- 
trated in a shell of width ~ 0.1i? m ; x . 

This contrasts with the survival of positrons in a stan- 
dard SN, where a substantial number of positrons can be 
created through the decay chain 56 Ni— > 56 Co^ 56 Fe+e + , 
which has a decay time of r ~ 111 days and occurs for 
18% of the 56 Ni atoms. [In SN1987A, observations imply 
that ~ 2.6 x 10 53 e + were produced through this chan- 
nel (Suntzeff & Bouchet 1990).] The SN positrons are 
produced inside the cool, dense SN ejecta when r ann /< <; 
0.03 (depending on the uncertain ionization and temper- 
ature structure of the ejecta). Therefore, most of the SN 
positrons are annihilated within a few years of the SN 
explosion. Milne, The, & Leising (1999) have confirmed 
numerically that only a small fraction of the positrons can 
escape the initial stages of the remnant evolution in Type 
la SNe. 

After passing through the external shock most positrons 
will rapidly thermalize with the ambient medium while 
a small fraction (~ 10 -3 ) will be re-accelerated by the 
(collisionless) external shock (Gieseler, Jones, & Kang 
2000). If positrons are produced ahead of the shock, then 
they will have a characteristic Lorentz factor 7+ ~ 30 
(Meszaros et al. 2001). Regardless of their initial state, 
the positrons cool rapidly through Coulomb collisions and 
adiabatic decompression as the remnant expands. While 
they are relativistic, adiabatic cooling leads to energy loss 
on the dynamical timescale of the remnant, r a( j ~ R s /v s , 
where R s and v s are the radius and expansion speed of 
the (assumed spherical) remnant. Once the positrons be- 
come non-relativistic, their Coulomb cooling time is short, 
Tboul ~ 3 x 10 4 (w+/c)rt ( j" 1 yr, where i>+ is the positron ve- 
locity. Using the analytic expressions for the cooling rates 
in Furlanetto & Loeb (2002) supplemented by expansion 
cooling, we have confirmed that the cooling timescales for 
7 + <J 100 are much shorter than the duration of the Sedov 
phase. We therefore simply assume that all positrons be- 
gin in thermal equilibrium with the shocked ambient gas. 

Once the remnant isotropizes, it follows the Sedov 
(1959) solution until radiative cooling becomes important. 
Although the remnant is initially asymmetric because of 
the jet geometry, Ayal & Piran (2001) have shown nu- 
merically that even before isotropization, the total vol- 
ume of the remnant scales with time similarly to a Se- 
dov self-similar blast wave. For this reason, and be- 
cause i; so is much smaller than the duration of the adi- 
abatic phase, we treat the remnant boundary as a Se- 
dov blast wave throughout both of these phases. Then 

R s = 0.314(i?5i/7io) 1 ^ 5 iyr 5 pc, where t yr is the remnant 
age in years. The adiabatic index of the ambient medium 
is taken to be T — 5/3. 

To calculate N+(t), we assume that that the positrons 



are perfectly mixed throughout the interior of the remnant, 
so that the e + density n + (r) ex n e (r), where n e is the free 
electron density. We expect such mixing to occur on a 
timescale ~ t- lso or faster because of instabilities on the 
interface between the SN ejecta and the ambient medium 
(Chevalier, Blondin, & Emmering(1992); Jun, Jones, & 
Norman 1996). The annihilation rate is then 

f R. 

N + = - I n e (r)n + {r)a[T(r)]tt-Kr 2 dr (1) 
Jo 

where the radial distributions of n e and T are given by 
Sedov (1959). Here a — ctf a + ap s is the annihilation 
rate coefficient with contributions from two channels. The 
first component, af a , involves direct annihilation with free 
electrons. In the non-relativistic regime (T <; 10 8 K), we 
use the rate coefficient given by Gould (1989) which in- 
cludes corrections due to Coulomb focusing. Above this 
temperature we use the fitting formula given by Svensson 
(1982). The second coefficient ap s describes radiative re- 
combination to positronium (Ps) and is given by Gould 
(1989). Although the corresponding expression is strictly 
valid only in the non-relativistic regime, Ps formation is 
sub-dominant for T ann > 10 6 K and so the errors intro- 
duced by applying it to high temperatures are negligible. 
The Ps lifetime is Tp s < 1.33 x 10~ 7 sec (depending on the 
energy level to which recombination occurs), and so we as- 
sume that positrons annihilate immediately after forming 
Ps. 

Once the shocked gas cools radiatively to T ^ 10 4 K, 
a dense shell begins to form at the shock front. We 
use the cooling function of Cioffi, McKee, & Bertschinger 
(1988), for which the nominal shell formation time is 

t s { = 3.61 x 10 4 -Egi n 4 ^ 7 yr. The shell formation process 
is not fully understood, but two-dimensional simulations 
have shown it to be both unstable and violent (Blondin 
et al. 1998). These simulations also find that the dense 
shell does not fully form until ~ 1.5t${. We therefore show 
results up to this time with the caveat that once cooling 
becomes important at t s t they may no longer be accurate. 
After shell formation we expect the positrons to annihilate 
rapidly. Not only does the density of the shell increase 
dramatically, but also charge exchange with HI comes to 
dominate as the temperature falls below ~ 10 5 K (Bus- 
sard, Ramaty, & Drachman 1979; see §^). Thus, we do 
not expect the annihilation signal to persist after the Se- 
dov phase ends. 

Even before shell formation, the loss of energy to cooling 
radiation decelerates the remnant. We use the analytic 
approximation of Cioffi et al. (1988) to describe the blast 
wave size and velocity in this phase (for t > t s {/e, where e 
is Euler's constant). As there is no analytic understanding 
of the interior structure of the remnant at this stage, we 
continue to use profiles from the Sedov (1959) solution, 
normalized to the total size and velocity of Cioffi et al. 
(1988). Our results are only weakly dependent on the 
details of these solutions. 

3. RESULTS 

Figure [j] shows the expected flux of annihilation pho- 
tons, Fj, as a function of time for blast waves traveling 
through media with densities tiq = 10 3 , 10 2 , 10 1 , 1, and 
0.1 cm~ 3 (from top to bottom). For each density we carry 



3 



the integration through 1.5t s f. The solid lines show the 
photon flux from direct annihilation while the dashed lines 
show that from Ps formation. We scale our results to d±o, 
the remnant distance in units of 10 kpc. The photon pro- 
duction rate is nearly proportional to the ambient den- 
sity because at a fixed temperature r ann oc n^ 1 . We find 
that the average temperature of annihilating positrons is 
T allI1 « 1.25T S , where T s is the postshock temperature. 
The annihilation weighted density is n ann ~ 2.7no. Both 
of these formulae are valid to within ~ 5% for all times 
and densities. Thus our model finds that most of the anni- 
hilations occur in a shell near the shock front. Throughout 
most of the remnant lifetime, 5 x 10 5 K <; T ann <; 10 7 K. 
Thus, direct e + e~ annihilation dominates, with Ps forma- 
tion becoming important only at late times when the shock 
velocity has fallen substantially. Note that F 7 is approxi- 
mately proportional to E51, but the explosion energy also 
affects t s {. 
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Fig. 1. — Flux of annihilation photons as a function of time for 
remnants embedded in a medium of density no = 10 3 , 10 2 , 10, 1, and 
0.1 cm" 3 , from top to bottom. Solid curves show the photon flux 
from free annihilation, while dashed curves show the flux from Ps 
decay. All results are scaled to a remnant distance dio = (d/10 kpc) 
and assume E51 = 1. 

We next calculate the emission spectrum of the rem- 
nant. The spectrum produced by direct annihilation when 
ksT e <C m e c 2 is given by Svensson, Larsson, & Poutanen 
(1996). We ignore Coulomb corrections to the annihila- 
tion rate, which is a good assumption for the temperatures 
of interest here (Gould 1989). The decay of Ps is more 
complicated. One-quarter of the recombinations produce 
para-Ps, which annihilates into two photons with the spec- 
tral distribution given by Gould (1989). Three-quarters of 
the recombinations produce ortho-Ps, which must decay 
into three photons in order to conserve angular momen- 
tum. The spectrum in this case is a continuum, with a 
form in the Ps rest frame given by Ore & Powell (1949). 
We include the broadening due to the thermal motions of 
the plasma particles by transforming to the observer frame 
(Bussard et al. 1979) and usingthe relativistic Maxwell- 
Boltzmann velocity distribution]]] In calculating the spec- 

1 Note, however, that our expression for «p s does not include rel- 
ativistic corrections, so our continuum level may be an overestimate 



tra, we assume for simplicity that all annihilations take 
place at a single temperature T ann . We also neglect the 
broadening introduced by the expansion speed of the rem- 
nant, because the ratio between the thermal e + speed and 
the shock speed is v^/vg ~ (m p /2m e ) 1 / 2 3> 1, where m e 
and m p are the electron and proton masses. 

Snapshots of the annihilation spectrum at t/t s f = 0.3, 0.6, 0.9, 1.2, 
and 1.5 are shown in Figure ^, for media with no = 1 cm" 3 
(top panel) and uq — 10 3 cm -3 (bottom panel). As shown 
in Figure [j], the number of free annihilations per sec- 
ond varies relatively slowly with remnant age, increasing 
primarily because of the decreasing T ann . The line also 
narrows with time: the full width at half maximum is 
~ 30(T ann /10 7 K) 1 / 2 keV. The line is narrower for lower 
density media at constant t/t s { because v s and hence T ann 
are smaller. 
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Fig. 2. — Spectral fluxes as a function of photon energy at 
various times. Results are shown for no = 1 cm -3 (top panel) and 
no = 10 3 cm" 3 (bottom panel) at times t/t B f = 0.3 (dotted curves), 
0.6 (dot-dashed curves), 0.9 (short-dashed curves), 1.2 (long-dashed 
curves), and 1.5 (solid curves). All results are scaled to a remnant 
distance dio = (d/10 kpc) and assume Est = 1. The asymmetric 
line broadening is a property of the annihilation spectrum given by 
Svensson et al. (1996). 

4. DISCUSSION 

The BATSE data set implies a GRB rate per comoving 
volume of hi so ~ 0.5 Gpc~ 3 yr" 1 , assuming isotropic emis- 
sion (Schmidt 2001). The actual event rate is estimated 
to be larger by an average beaming factor, (/ b _1 ) ~ 520 
(Frail et al. 2001). Most models associate GRBs with 
compact stellar remnants; under the assumption that 
the GRB event rate is proportional to the star forma- 
tion rate we can convert the above rate density to an 
event rate per unit mass of star formation (e.g., Por- 
ciani & Madau 2001). The expected period between 
GRB events in the Milky Way is then tmw ~ 3 x 

10 4 yr (n iso /0.5)- 1 ((/ ( T 1 )/520)- 1 (i? /0.007)i 1 , M 1 W' where 
i?o and Rmw are the present day average star formation 
rates in the universe and in the Milky Way in units of 
M Q Mpc -3 yr -1 and M yr -1 , respectively, and ni SO has 
units of Gpc -3 yr -1 . 

at very early times when T ann > 10 s K. 
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Together with Figure [l], this suggests that the positron 
annihilation line should be detectable from at least one 
Galactic GRB remnant at any time provided that the cir- 
cumburst medium has no <J 1 cm -3 . This conclusion re- 
lies on three fortunate coincidences. First and most im- 
portantly, the shell formation time satisfies coincidentally 
t s f ~ Tmw for no = 1 cm~ 3 . Although some fraction of 
GRBs may occur in dense molecular clouds, recent mod- 
eling of jetted afterglows indicates that the ambient den- 
sity is in the range no ~ 10 _3 -10 15 cm -3 (Berger et al. 
2000, 2001; Harrison et al. 2001; Panaitescu & Kumar 
2001a). Second, in contrast to e + created by SNe, GRB 
positrons annihilate on a timescale r ann > t s f. Third, 

Tknn > 5 x 10 5 K, so that a clear line from the direct 
annihilation channel is visible throughout the source life- 
time. 

Two backgrounds could contaminate the annihilation 
line signal. First, relativistic protons accelerated by the 
collisionless shock can create positrons through the de- 
cay of 7r + produced in collisons with thermal protons. 
Assuming that a fraction £ p of the GRB energy is in- 
vested in the acceleration of protons (Blandford & Eich- 
ler 1987) and using the approximate interaction cross- 
sections of Mannheim & Schlickeiser (1994), we find that 
the number of positrons produced before shell forma- 
tion is ~ lQ 48 Eg[ /14 no /7 (£ P /0.1) (see Furlanetto & 
Loeb 2002). Unless £-2.3 ^ 10~ 6 , the positron popula- 
tion from this process is negligible compared to that pro- 
duced in the GRB. A second background is due to inverse- 
Compton (IC) emission by the electrons accelerated at 
the shock. The recent detection of TeV 7-ray emission 
from SN remnants (Tanimori et al. 1998; Muraishi et 
al. 2000) coupled with the well-observed high-energy syn- 
chrotron emission from the same remnants (Koyama et 
al. 1995, 1997) provides strong evidence for the existence 
of this population of electrons and calibrates the IC sig- 
nal. If we extrapolate the observed TeV 7-ray flux to the 
energy of the annihilation line using the spectral index 
of the radio synchrotron emission (which should also be 
the index of the IC spectrum; Reynolds 2001), we find a 
flux F 7) io < 10- 12 di" 2 (e 7 /511 keV)" 15 cm" 2 s" 1 keV" 1 . 
Thus, the IC background should be negligible compared 
to the annihilation signal. 

The prospects for detection of young GRB remnants 
in the Milky Way galaxy are good, provided that GRB 
sources are not confined to dense molecular clouds. As 
shown in Figure [l], a GRB remnant produces a char- 
acteristic annihilation line photon flux of F y ~ 3 x 
10~ 5 £-2.3no/dio cm" 2 s -1 . The INTEGRAL satellite^ 
expected to be launched in October 2002, will have spec- 
tral capabilities in the energy range of interest. For a 
10 6 sec observation, its SPI and IBIS instruments have 
3cr line sensitivities of ~ 5 x 10~ 6 cm -2 s _1 and ~ 2 x 
10~ 5 cm -2 s _1 , respectively. EXIST, a proposed all-sky 
hard X-ray survey mission, has an expected 5<r line sensi- 
tivity of ~ 5 x 10~ 6 cm~ 2 s _1 in the relevant energy range 
(assuming an integration time of 10 7 sec, the mean expo- 
sure time planned for any point on the sky). These instru- 
ments will also have the spectral resolution to identify the 
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line. The IBIS instrument on INTEGRAL and the EXIST 
survey, with angular resolutions of 12' and 5' respectively, 
will also be able to map nearby GRB remnants, for which 
the angular size at t s { is ~ lAE^p n 3 ^ 7 d^Q arcmin. 

If GRBs occur primarily in dense environments, t a t de- 
creases and the probability for observing a remnant de- 
clines. Figure [l] shows that the average flux in the an- 
nihilation line from a nearby extragalactic source (d ;> 
1 Mpc) during the Sedov phase will likely remain below 
the sensitivity limits of the above instruments. How- 
ever, once shell formation begins, the positrons will an- 
nihilate through charge exchange with recombining hy- 
drogen on a timescale r cx ~ 10/(/Hi^ s h) y r , where /hi 
is the neutral fraction (Bussard et al. 1979). If the 
shell is supported primarily by magnetic and cosmic ray 
pressure, Cox et al. (1999) find that the shell density 

at t s { is rt s h ~ IOQBq 1 !!^ 3 ^ 14 El{ 14 cm~ 3 , where B is 
the transverse component of the ambient magnetic field 
in units of /zG. The maximal annihilation flux is then 
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where dM pc is the distance to the GRB in Mpc. Because 
the charge exchange process produces Ps, 75% of the an- 
nihilation energy is released in the continuum, but this 
suggests that the remnants will brighten significantly dur- 
ing the brief phase of shell formation and may be visible 
out to the Virgo cluster distance, d ~ 20 Mpc. Given the 
number of galaxies in Virgo, we would therefore expect 
to find at least one Virgo remnant in this stage provided 
that fmn s h <, 1 cm~ 3 . The details of the luminosity and 
spectrum will depend critically on modeling the process 
of shell formation and cooling. Such modeling stretches 
the limits of current simulations (e.g., Blondin et al. 1998; 
Cox et al. 1999) but may become feasible in the future. 
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was supported in part by NASA grants NAG 5-7039, 5- 
7768, and NSF grants AST-9900877, AST-0071019 for AL. 
SRF acknowledges the support of an NSF graduate fellow- 
ship. 



REFERENCES 

Ayal, S. & Piran, T. 2001, ApJ, 555, 23 

Berger, E. et al. 2000, ApJ, 545, 56 

Berger, E. et al. 2001, ApJ, 556, 556 

Blandford, R. D. & Eichler, D. 1987, Phys. Rep., 154, 1 

Blondin, J. M., Wright, E. B., Borkowski, K. J., & Reynolds, S. P. 

1998, ApJ, 500, 342 
Bloom, J. S. et al. 1999, Nature, 401, 453 

Bussard, R. W., Ramaty, R., & Drachman, R. J. 1979, ApJ, 228, 
928 

Cavallo, G. & Rees, M. J. 1978, MNRAS, 183, 359 
Chevalier, R. A., Blondin, J. M., & Emmering, R. T. 1992, ApJ, 392, 
118 

Cioffi, D. F., McKee, C. F., & Bertschinger, E. 1988, ApJ, 334, 252 
Cox, D. P. et al. 1999, ApJ, 524, 179 

Dermer, C. D. & Bottcher, M. 2000, ApT, LI55 

Djorgovski, S. C, et al. 2001, preprint (a.stro-ph/0107535) 

Frail, D. A. et al. 2001, ApJ, 562, L55 

Freedman, D. L. & Waxman, E. 2001, ApJ, 547, j932 . 

Furlanetto, S. R. & Loeb, A. 2002, ApJ, in press ( astro-ph/0201313 ) 
Gieseler, U. D. J., Jones, T. W., & Kang, H. 2000, A&A, 364, 911 
Gould, R. J. 1989, ApJ, 344, 232 
Harrison, F. A. et al. 2001, ApJ, 559, 123 

Jun, B., Jones, T. W., & Norman, M. L. 1996, ApJ, 468, L59 



Koyama, K. et al. 1995, Nature, 378, 255 
Koyama, K. et al. 1997, PASJ, 49, L7 



5 



Kulkarni, S. R., et al. 2000, Proc. SPIE, 4005, 9 

Loeb, A. & Perna, R. 1998, ApJ, 503, L35 

Mannheim, K. & Schlickeiser, R. 1994, A&A, 286, 983 

Meszaros, P., Ramirez-Ruiz, E., & Rees, M. J. 2001, ApJ, 554, 660 

Milne, P. A., The, L.-S., & Leising, M. D. 1999, ApJS, 124, 503 

Muraishi, H. et al. 2000, A&A, 354, L57 

Ore, A. & Powell, J. L. 1949, Phys Rev, 75, 1696 

Paczynski, B. 2001, Acta Astronomica, 51, 1 

Panaitescu, A. & Kumar, P. 2001a, ApJ, 554, 667 

. 2001b, ApJ, 560, L49 

Perna, R., Raymond, J., & Loeb, A. 2000, ApJ, 533, 658 
Piran, T. 2000, Phys. Rep., 333, 529 

Piran, T. & Ayal, S. 2002, in Disks of Galaxies: Kinematics, 

Dynamics and Interactions, eds. Athaaasstmla, E & Bosma, A., 

to appear in ASP Conf. Series (a.stro-ph/0203003) 

Piran, T., Kumar, P., Panaitescu, A., & Piro, L. 2001, ApJ, 560, 
L167 

Porciani, C. & Madau, P. 2001, ApJ, 548, 522 

Reichart, D. E. 2001, ApJ, 554, 643 

Reynolds, S. P. 2001, Space Science Reviews, 99, 177 

Rhoads, J. E. 1997, ApJ, 487, LI 

Schmidt, M. 2001, ApJ, 552, 36 

Sedov, L. I. 1959, Similarity and Dimensional Methods in Mechanics 

(New York: Academic Press) 
Shemi, A. & Piran, T. 1990, ApJ, 365, L55 
Suntzeff, N. B. & Bouchet, P. 1990, AJ, 99, 650 
Svensson, R., Larsson, S., & Poutanen, J. 1996, A&AS, 120, C587 
Svensson, R. 1982, ApJ, 258, 32 
Tanimori, T. et al. 1998, ApJ, 497, L25 
Thompson, C. & Madau, P. 2000, ApJ, 538, 105 
Truelove, J. K. & McKee, C. F. 1999, ApJS, 120, 299 
van Paradijs, J., Kouveliotou, C, & Wijers, R. A. M. J. 2000, 
ARA&A, 38, 379 



Woods, E., & Loeb, A. 1999, astro-ph/9907110 



